Maturation and Conformational Switching of a De Novo Designed Phase-Separating Polypeptide.

Cellular compartments formed by biomolecular condensation are widespread features of cell biology. These organelle-like assemblies compartmentalize macromolecules dynamically within the crowded intracellular environment. However, the intermolecular interactions that produce condensed droplets may also create arrested states and potentially pathological assemblies such as fibers, aggregates, and gels through droplet maturation. Protein liquid-liquid phase separation is a metastable process, so maturation may be an intrinsic property of phase-separating proteins, where nucleation of different phases or states arises in supersaturated condensates. Here, we describe the formation of both phase-separated droplets and proteinaceous fibers driven by a de novo designed polypeptide. We characterize the formation of supramolecular fibers in vitro and in bacterial cells. We show that client proteins can be targeted to the fibers in cells using a droplet-forming construct. Finally, we explore the interplay between phase separation and fiber formation of the de novo polypeptide, showing that the droplets mature with a post-translational switch to largely ß conformations, analogous to models of pathological phase separation.

Modeling Zinc Complexes Using Neural Networks.

Understanding the energetic landscapes of large molecules is necessary for the study of chemical and biological systems. Recently, deep learning has greatly accelerated the development of models based on quantum chemistry, making it possible to build potential energy surfaces and explore chemical space. However, most of this work has focused on organic molecules due to the simplicity of their electronic structures as well as the availability of data sets. In this work, we build a deep learning architecture to model the energetics of zinc organometallic complexes. To achieve this, we have compiled a configurationally and conformationally diverse data set of zinc complexes using metadynamics to overcome the limitations of traditional sampling methods. In terms of the neural network potentials, our results indicate that for zinc complexes, partial charges play an important role in modeling the long-range interactions with a neural network. Our developed model outperforms semiempirical methods in predicting the relative energy of zinc conformers, yielding a mean absolute error (MAE) of 1.32 kcal/mol with reference to the double-hybrid PWPB95 method.

Computational insights into the conformational transition of STING: Mechanistic, energetic considerations, and the influence of crucial mutations.

STING (stimulator of interferon genes) is a crucial protein in the innate immune system's response to viral and bacterial infections. In this study, we investigated the mechanistic and energetic mechanism of the conformational transition process of STING activated by cGAMP binding. We found that the STING connector region undergoes an energetically unfavorable rotation during this process, which is compensated by the favorable interaction between cGAMP and the STING ligand binding domain. We further studied several disease-causing mutations and found that the V155 M mutation facilitates a smoother transition in the STING connector region. However, the V147L mutation exhibits unfavorable conformational transition energy, suggesting it may hinder STING activation pathway that relies on connector region rotation. Despite being labeled as hyperactive, the widespread prevalence of V147L/V147I mutations across species implies a neutral character, indicating complexity in its role. Overall, our analysis deepens the understanding of STING activation within the connector region, and targeting this region with compounds may provide an alternative approach to interfering with STING's function.

Molecular dynamics simulation of the effect of temperature on the conformation of ubiquitin protein.

CONTENT: Ubiquitin, a ubiquitous small protein found in all living organisms, is crucial for tagging proteins earmarked for degradation and holds pivotal importance in biomedicine. Protein functionality is intricately linked to its structure. To comprehend the impact of diverse temperatures on ubiquitin protein structure, our study delved into the energy landscape, hydrogen bonding, and overall structural stability of ubiquitin protein at varying temperatures. Through meticulous analysis of root mean square deviation and root mean square fluctuation, we validated the robustness of the simulation conditions employed. Within our simulated system, the bonding energy and electrostatic potential energy exhibited linear augmentation, while the van der Waals energy demonstrated a linear decline. Additionally, our findings highlighted that the α-Helix secondary structure of the ubiquitin protein gradually transitions toward helix destabilization under high-temperature conditions. The secondary structure of ubiquitin protein experiences distinct changes under varying temperatures. The outcomes of our molecular simulations offer a theoretical framework that enhances our comprehension of how temperature impacts the structural stability of ubiquitin protein. These insights contribute not only to a deeper understanding of iniquity's behavior but also hold broader implications in the realm of biomedicine and beyond. METHODS: All the MD simulations were performed using the GROMACS software with GROMOS96 force field and SPC for water. The ubiquitin protein was put in the center of a cubic box with a length of 8 nm, a setting that allowed > 0.8 nm in the minimal distance between the protein surface and the box wall. To remove the possible coordinate collision of the configurations, in the beginning, the steepest descent method was used until the maximum force between atoms was under 100 kJ/mol·nm with a 0.01 nm step size. Minimization was followed by 30 ps of position-restrained MD simulation. The protein was restrained to its initial position, and the solvent was freely equilibrated. The product phase was obtained with the whole system simulated for 10 ns without any restraint using an integral time step of 1 fs with different temperatures. The cutoff for short-range electronic interaction was set to 1.5 nm. The long-range interactions were treated with a particle-mesh Ewald (PME) method with a grid width of 1.2 nm.

Conformational Effects on the Absorption Spectra of Phytochromes.

Bacteriophytochrome is a photoreceptor protein that contains the biliverdin (BV) chromophore as its active component. The spectra of BV upon mutation remain remarkably unchanged, as far as spectral positions are concerned. This points toward the minimal effect of electrostatic effects on the electronic structure of the chromophore. However, the relative intensities of the Q and Soret bands of the chromophore change dramatically upon mutation. In this work, we delve into the molecular origin of this unusual intensity modulation. Using extensive classical MD and QM/MM calculations, we show that due to mutation, the conformational population of the chromophore changes significantly. The noncovalent interactions, especially the stacking interactions, lead to extra stabilization of the cyclic form in the D207H mutated species as opposed to the open form in the wild-type BV. Thus, unlike the commonly observed direct electrostatic effect on the spectral shift, in the case of BV the difference observed is in varying intensities, and this in turn is driven by a conformational shift due to enhanced stacking interaction.

Photophysics of a nucleic acid-protein crosslinking model strongly depends on solvation dynamics: an experimental and theoretical study.

We present a combined experimental and theoretical study of the photophysics of 5-benzyluracil (5BU) in methanol, which is a model system for interactions between nucleic acids and proteins. A molecular dynamics study of 5BU in solution through efficient DFT-based hybrid ab initio potentials revealed a remarkable conformational flexibility - allowing the population of two main conformers - as well as specific solute-solvent interactions, which both appear as relevant factors for the observed 5BU optical absorption properties. The simulated absorption spectrum, calculated on such an ensemble, enabled a molecular interpretation of the experimental UV-Vis lowest energy band, which is also involved in the induced photo-reactivity upon irradiation. In particular, the first two excited states (mainly involving the uracil moiety) both contribute to the 5BU lowest energy absorption. Moreover, as a key finding, the nature and brightness of such electronic transitions are strongly influenced by 5BU conformation and the microsolvation of its heteroatoms.

Tertiary structure and conformational dynamics of the anti-amyloidogenic chaperone DNAJB6b at atomistic resolution.

DNAJB6b is a molecular chaperone of the heat shock protein network, shown to play a crucial role in preventing aggregation of several disease-related intrinsically disordered proteins. Using homology modeling and microsecond-long all-atom molecular dynamics (MD) simulations, we show that monomeric DNAJB6b is a transiently interconverting protein cycling between three states: a closed state, an open state (both abundant), and a less abundant extended state. Interestingly, the reported regulatory autoinhibitory anchor between helix V in the G/F1 region and helices II/III of the J-domain, which obstructs the access of Hsp70 to the J-domain remains present in all three states. This possibly suggests a mechanistically intriguing regulation in which DNAJB6b only becomes exposed when loaded with substrates that require Hsp70 processing. Our MD results of DNAJB6b carrying mutations in the G/F1 region that are linked to limb-girdle muscular dystrophy type D1 (LGMDD1) show that this G/F1 region becomes highly dynamic, pointing towards a spontaneous release of the autoinhibitory helix V from helices II/III. This would increase the probability of non-functional Hsp70 interactions to DNAJB6b without substrates. Our cellular data indeed confirm that non-substrate loaded LGMDD1 mutants have aberrant interactions with Hsp70.

Molecular dynamics in multidimensional space explains how mutations affect the association path of neomycin to a riboswitch.

Riboswitches are messenger RNA (mRNA) fragments binding specific small molecules to regulate gene expression. A synthetic N1 riboswitch, inserted into yeast mRNA controls the translation of a reporter gene in response to neomycin. However, its regulatory activity is sensitive to single-point RNA mutations, even those distant from the neomycin binding site. While the association paths of neomycin to N1 and its variants remain unknown, recent fluorescence kinetic experiments indicate a two-step process driven by conformational selection. This raises the question of which step is affected by mutations. To address this, we performed all-atom two-dimensional replica-exchange molecular dynamics simulations for N1 and U14C, U14C[Formula: see text], U15A, and A17G mutants, ensuring extensive conformational sampling of both RNA and neomycin. The obtained neomycin association and binding paths, along with multidimensional free-energy profiles, revealed a two-step binding mechanism, consisting of conformational selection and induced fit. Neomycin binds to a preformed N1 conformation upon identifying a stable upper stem and U-turn motif in the riboswitch hairpin. However, the positioning of neomycin in the binding site occurs at different RNA-neomycin distances for each mutant, which may explain their different regulatory activities. The subsequent induced fit arises from the interactions of the neomycin's N3 amino group with RNA, causing the G9 backbone to rearrange. In the A17G mutant, the critical C6-A17/G17 stacking forms at a closer RNA-neomycin distance compared to N1. These findings together with estimated binding free energies coincide with experiments and elucidate why the A17G mutation decreases and U15A enhances N1 activity in response to neomycin.

Using Three-Dimensional Information to Predict and Interpret the Facial Selectivities of Nucleophilic Additions to Cyclic Ketones.

In this study, we devised a new method to predict facial selectivity by quantifying steric and orbital factors for the nucleophile approaching both π-plane faces. Using this method, we quantified the total electron density and frontier orbital distributions of 163 cyclic ketones with various structures and quantitatively explained the surface selectivity of 323 reactions with eight nucleophiles (BH3, LiAlH4, NaBH4, LiAl(OMe)3H, MeLi, MeMgI, PhLi, and PnMgI). Importance analysis showed a large orbital effect for BH3, LiAlH4, and NaBH4 and the dominance of the steric effect for LiAl(OMe)3H, MeLi, MeMgI, PhLi, and PhMgI. Our method analyzes three-dimensional features based on Gaussian cube files, which can be easily obtained using mainstream computational chemistry software packages, and this approach should prove useful for predicting the rates and facial selectivity of other reactions.

Transient interdomain interactions in free USP14 shape its conformational ensemble.

The deubiquitinase (DUB) ubiquitin-specific protease 14 (USP14) is a dual domain protein that plays a regulatory role in proteasomal degradation and has been identified as a promising therapeutic target. USP14 comprises a conserved USP domain and a ubiquitin-like (Ubl) domain separated by a 25-residue linker. The enzyme activity of USP14 is autoinhibited in solution, but is enhanced when bound to the proteasome, where the Ubl and USP domains of USP14 bind to the Rpn1 and Rpt1/Rpt2 units, respectively. No structure of full-length USP14 in the absence of proteasome has yet been presented, however, earlier work has described how transient interactions between Ubl and USP domains in USP4 and USP7 regulate DUB activity. To better understand the roles of the Ubl and USP domains in USP14, we studied the Ubl domain alone and in full-length USP14 by nuclear magnetic resonance spectroscopy and used small angle x-ray scattering and molecular modeling to visualize the entire USP14 protein ensemble. Jointly, our results show how transient interdomain interactions between the Ubl and USP domains of USP14 predispose its conformational ensemble for proteasome binding, which may have functional implications for proteasome regulation and may be exploited in the design of future USP14 inhibitors.

Expansive discovery of chemically diverse structured macrocyclic oligoamides.

Small macrocycles with four or fewer amino acids are among the most potent natural products known, but there is currently no way to systematically generate such compounds. We describe a computational method for identifying ordered macrocycles composed of alpha, beta, gamma, and 17 other amino acid backbone chemistries, which we used to predict 14.9 million closed cycles composed of >42,000 monomer combinations. We chemically synthesized 18 macrocycles predicted to adopt single low-energy states and determined their x-ray or nuclear magnetic resonance structures; 15 of these were very close to the design models. We illustrate the therapeutic potential of these macrocycle designs by developing selective inhibitors of three protein targets of current interest. By opening up a vast space of readily synthesizable drug-like macrocycles, our results should considerably enhance structure-based drug design.

Synthesis, crystal structure and in-silico evaluation of arylsulfonamide Schiff bases for potential activity against colon cancer.

This report presents a comprehensive investigation into the synthesis and characterization of Schiff base compounds derived from benzenesulfonamide. The synthesis process, involved the reaction between N-cycloamino-2-sulfanilamide and various substituted o-salicylaldehydes, resulted in a set of compounds that were subjected to rigorous characterization using advanced spectral techniques, including 1H NMR, 13C NMR and FT-IR spectroscopy, and single-crystal X-ray diffraction. Furthermore, an in-depth assessment of the synthesized compounds was conducted through Absorption, Distribution, Metabolism, Excretion and Toxicity (ADMET) analysis, in conjunction with docking studies, to elucidate their pharmacokinetic profiles and potential. Impressively, the ADMET analysis showcased encouraging drug-likeness properties of the newly synthesized Schiff bases. These computational findings were substantiated by molecular properties derived from density functional theory (DFT) calculations using the B3LYP/6-31G* method within the Jaguar Module of Schrödinger 2023-2 from Maestro (Schrodinger LLC, New York, USA). The exploration of frontier molecular orbitals (HOMO and LUMO) enabled the computation of global reactivity descriptors (GRDs), encompassing charge separation (Egap) and global softness (S). Notably, within this analysis, one Schiff base, namely, 4-bromo-2-{N-[2-(pyrrolidine-1-sulfonyl)phenyl]carboximidoyl}phenol, 20, emerged with the smallest charge separation (ΔEgap = 3.5780 eV), signifying heightened potential for biological properties. Conversely, 4-bromo-2-{N-[2-(piperidine-1-sulfonyl)phenyl]carboximidoyl}phenol, 17, exhibited the largest charge separation (ΔEgap = 4.9242 eV), implying a relatively lower propensity for biological activity. Moreover, the synthesized Schiff bases displayed remarkeable inhibition of tankyrase poly(ADP-ribose) polymerase enzymes, integral in colon cancer, surpassing the efficacy of a standard drug used for the same purpose. Additionally, their bioavailability scores aligned closely with established medications such as trifluridine and 5-fluorouracil. The exploration of molecular electrostatic potential through colour mapping delved into the electronic behaviour and reactivity tendencies intrinsic to this diverse range of molecules.

Sequence-dependent material properties of biomolecular condensates and their relation to dilute phase conformations.

Material properties of phase-separated biomolecular condensates, enriched with disordered proteins, dictate many cellular functions. Contrary to the progress made in understanding the sequence-dependent phase separation of proteins, little is known about the sequence determinants of condensate material properties. Using the hydropathy scale and Martini models, we computationally decipher these relationships for charge-rich disordered protein condensates. Our computations yield dynamical, rheological, and interfacial properties of condensates that are quantitatively comparable with experimentally characterized condensates. Interestingly, we find that the material properties of model and natural proteins respond similarly to charge segregation, despite different sequence compositions. Molecular interactions within the condensates closely resemble those within the single-chain ensembles. Consequently, the material properties strongly correlate with molecular contact dynamics and single-chain structural properties. We demonstrate the potential to harness the sequence characteristics of disordered proteins for predicting and engineering the material properties of functional condensates, with insights from the dilute phase properties.

Role of Protonation States in the Stability of Molecular Dynamics Simulations of High-Resolution Membrane Protein Structures.

Classical molecular dynamics (MD) simulations provide unmatched spatial and time resolution of protein structure and function. However, the accuracy of MD simulations often depends on the quality of force field parameters and the time scale of sampling. Another limitation of conventional MD simulations is that the protonation states of titratable amino acid residues remain fixed during simulations, even though protonation state changes coupled to conformational dynamics are central to protein function. Due to the uncertainty in selecting protonation states, classical MD simulations are sometimes performed with all amino acids modeled in their standard charged states at pH 7. Here, we performed and analyzed classical MD simulations on high-resolution cryo-EM structures of two large membrane proteins that transfer protons by catalyzing protonation/deprotonation reactions. In simulations performed with titratable amino acids modeled in their standard protonation (charged) states, the structure diverges far from its starting conformation. In comparison, MD simulations performed with predetermined protonation states of amino acid residues reproduce the structural conformation, protein hydration, and protein-water and protein-protein interactions of the structure much better. The results support the notion that it is crucial to perform basic protonation state calculations, especially on structures where protonation changes play an important functional role, prior to the launch of any conventional MD simulations. Furthermore, the combined approach of fast protonation state prediction and MD simulations can provide valuable information about the charge states of amino acids in the cryo-EM sample. Even though accurate prediction of protonation states in proteinaceous environments currently remains a challenge, we introduce an approach of combining pKa prediction with cryo-EM density map analysis that helps in improving not only the protonation state predictions but also the atomic modeling of density data.

Conformational Heterogeneity of ß-Barrel Membrane Proteins Observed In Situ Using Orthogonal Spin Labels and Pulsed ESR Spectroscopy.

Outer membrane proteins (OMPs) of Gram-negative bacteria are involved in many essential functions of the cell. They are tightly packed in the outer membrane, which is an asymmetric lipid bilayer. Electron spin resonance (ESR) spectroscopic techniques combined with site-directed spin labeling (SDSL) enable observation of structure and conformational dynamics of these proteins directly in their native environments. Here we depict a protocol for site-directed spin labeling of ß-barrel membrane proteins in isolated outer membranes and intact E. coli using nitroxide, triarylmethyl (trityl), and Gd3+-based spin tags. Furthermore, subsequent continuous wave (CW) and orthogonal pulsed electron-electron double resonance (PELDOR) measurements are described along with experimental setup at Q-band (34 GHz), the data analysis, and interpretation.

Transient Receptor Potential Melastatin 7 (TRPM7) Ion Channel Inhibitors: Preliminary SAR and Conformational Studies of Xenicane Diterpenoids from the Hawaiian Soft Coral Sarcothelia edmondsoni.

Waixenicin A, a xenicane diterpene from the octocoral Sarcothelia edmondsoni, is a selective, potent inhibitor of the TRPM7 ion channel. To study the structure-activity relationship (SAR) of waixenicin A, we isolated and assayed related diterpenes from S. edmondsoni. In addition to known waixenicins A (1) and B (2), we purified six xenicane diterpenes, 7S,8S-epoxywaixenicins A (3) and B (4), 12-deacetylwaixenicin A (5), waixenicin E (6), waixenicin F (7), and 20-acetoxyxeniafaraunol B (8). We elucidated the structures of 3-8 by NMR and MS analyses. Compounds 1, 2, 3, 4, and 6 inhibited TRPM7 activity in a cell-based assay, while 5, 7, and 8 were inactive. A preliminary SAR emerged showing that alterations to the nine-membered ring of 1 did not reduce activity, while the 12-acetoxy group, in combination with the dihydropyran, appears to be necessary for TRPM7 inhibition. The bioactive compounds are proposed to be latent electrophiles by formation of a conjugated oxocarbenium ion intermediate. Whole-cell patch-clamp experiments demonstrated that waixenicin A inhibition is irreversible, consistent with a covalent inhibitor, and showed nanomolar potency for waixenicin B (2). Conformational analysis (DFT) of 1, 3, 7, and 8 revealed insights into the conformation of waixenicin A and congeners and provided information regarding the stabilization of the proposed pharmacophore.

Molecular Insights into the Impact of Mutations on the Binding Affinity of Targeted Covalent Inhibitors of BTK.

Targeted covalent inhibitors (TCIs) have witnessed a significant resurgence in recent years, particularly in the kinase drug discovery field for treating diverse clinical indications. The inhibition of Bruton's tyrosine kinase (BTK) for treating B-cell cancers is a classic example where TCIs such as ibrutinib have had breakthroughs in targeted therapy. However, selectivity remains challenging, and the emergence of resistance mutations is a critical concern for clinical efficacy. Computational methods that can accurately predict the impact of mutations on inhibitor binding affinity could prove helpful in informing targeted approaches─providing insights into drug resistance mechanisms. In addition, such systems could help guide the systematic evaluation and impact of mutations in disease models for optimal experimental design. Here, we have employed in silico physics-based methods to understand the effects of mutations on the binding affinity and conformational dynamics of select TCIs of BTK. The TCIs studied include ibrutinib, acalabrutinib, and zanubrutinib─all of which are FDA-approved drugs for treating multiple forms of leukemia and lymphoma. Our results offer useful molecular insights into the structural determinants, thermodynamics, and conformational energies that impact ligand binding for this biological target of clinical relevance.

Genetically Fused Resilin-like Polypeptide-Coiled Coil Bundlemer Conjugates Exhibit Tunable Multistimuli-Responsiveness and Undergo Nanofibrillar Assembly.

Peptide-based materials are diverse candidates for self-assembly into modularly designed and stimuli-responsive nanostructures with precisely tunable compositions. Here, we genetically fused computationally designed coiled coil-forming peptides to the N- and C-termini of compositionally distinct multistimuli-responsive resilin-like polypeptides (RLPs) of various lengths. The successful expression of these hybrid polypeptides in bacterial hosts was confirmed through techniques such as gel electrophoresis, mass spectrometry, and amino acid analysis. Circular dichroism spectroscopy and ultraviolet-visible turbidimetry demonstrated that despite the fusion of disparate structural and responsive units, the coiled coils remained stable in the hybrid polypeptides, and the sequence-encoded differences in thermoresponsive phase separation of the RLPs were preserved. Cryogenic transmission electron microscopy and coarse-grained modeling showed that after thermal annealing in solution, the hybrid polypeptides adopted a closed loop conformation and assembled into nanofibrils capable of further hierarchically organizing into cluster structures and ribbon-like structures mediated by the self-association tendency of the RLPs.

Synthesis and ab initio conformational investigation of a series of model sulfated α-L-iduronopyranosides.

Due to the all-axial orientation of the OH-groups in the 1C4 chair conformation considered standard for L-hexapyranosides, including l-iduronopyranoside - a component of many biologically and medically significant sulfated glycans, these monosaccharides can be anticipated to display unusual conformations upon the introduction of bulky and charged substituents. Herein we describe the synthesis of a series of iduronopyranoside derivatives with varying sulfation patterns, which were studied computationally using the DLPNO-MP2 approach and by means of analyzing their chemical shifts to ascertain the effects sulfation has on the conformation of the iduronopyranoside ring.

General Synthesis of Conformationally Constrained Noncanonical Amino Acids with C(sp3)-Rich Benzene Bioisosteres.

Recent years have seen novel modalities emerge for the treatment of human diseases resulting in an increase in beyond rule of 5 (bRo5) chemical matter. As a result, synthetic innovations aiming to enable rapid access to complex bRo5 molecular entities have become increasingly valuable for medicinal chemists' toolkits. Herein, we report the general synthesis of a new class of noncanonical amino acids (ncAA) with a cyclopropyl backbone to achieve conformational constraint and bearing C(sp3)-rich benzene bioisosteres. We also demonstrate preliminary studies toward utilities of these ncAA as building blocks for medicinal chemistry research.